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I. INTRODUCTION 

In the past two decades quantum computation has at- 
tracted growing interest, encouraged by the development 
of tools for the manipulation of single quantum objects 
as well as by several remarkable theoretical findings [T]. 
Different systems have been proposed as candidates for 
quantum computing; they are based, for instance, on 
cavity-laser atoms, Bose-Einstein condensates, or NMR 
techniques (see e.g. At the same time, a num- 

ber of quantum algorithms have been designed and some 
have been shown to be even exponentially faster than 
their best classical counterparts [SI. In particular, quan- 
tum search algorithms, although able to achieve "only" a 
polynomial speedup, have been proved to be very promis- 
ing and of widespread use in quantum computation [5J[7]. 

One of the best known quantum search algorithms is 
attributed to Grover [5]. The algorithm can find a tar- 
get within an unsorted database made up of N items 
using 0{VN) queries. Due to its broad range of applica- 
tions and its ability to be effectively used as a subroutine 
in], Grover's algorithm has been thoroughly investi- 
gated and a number of different implementations have 
been proposed (see e.g. [H ISl [TUHTB] '). Recently, imple- 
mentations based on continuous-time [TTJ [T^ as well as 
on discrete-time quantum walks [HI UHl [H] have been 
introduced: Given that the application of random walks 
in classical algorithms provided significant advantages for 
approximations and optimization, one is strongly moti- 
vated to study quantum walks as algorithmic tools. At 
the current stage the aim is not only to achieve similar 
computational improvements, but also to understand the 
capabilities of quantum computations. 

Here we focus on the approach pioneered by Farhi 
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and Gutmann [TT] and further developed by Childs and 
Goldstone [T^] , based on continuous-time quantum walks 
(CTQWs). This implementation provides several impor- 
tant advantages: the algorithm does not need auxiliary 
storage space and it makes it possible to take into account 
the geometrical arrangement of the database (the latter 
being either a physical position space or an efficiently en- 
coded Hilbert space ;,14J). Indeed, while previous studies 
just considered the cases of the translationally invariant 
Z"* (named d-dimensional periodic lattices by Childs and 
Goldstone [T^j and called hypercubic lattices in statisti- 
cal physics, the term which we adopt here) and of com- 
plete graphs [21] , here we extend the investigations to the 
case of generic structures and analyze how geometrical 
parameters [e.g., (fractal) dimension or (average) coor- 
dination number; see later in this article] affect the dy- 
namics of the CTQW. In [12] it was shown that quantum 
searches based on CTQWs recover the optimal quadratic 
speedup on complete graphs and on high-dimensional hy- 
percubic lattices (with dimension rf > 4), while for low- 
dimensional (d < 4) lattices CTQWs can not outperform 
their classical counterpart. Hence, it may appear that 
d = 4 works as a "critical dimension" , separating highly 
performing structures from poorly performing ones. We 
will show, both analytically and numerically, that the 
dimension of the substrate is not sufficient for getting 
a sharp transition from the ground to the first excited 
state; we also take into account the success probability 
T^w.sit), that is, the probability of finding the quantum 
walker at the target site w at time t, given as initial 
state the equally weighted superposition s: An efficient 
quantum walk gives rise to a success probability close 
to 1 already at very small times t. In particular, we will 
take into account several kinds of structures: translation- 
ally invariant structures, such as d-dimensional hypercu- 
bic lattices, complete graphs, fractals with low fractal 
dimension like dual Sierpinski gaskets and T-fractals, hi- 
erarchical structures as Cayley trees, and structures with 
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(fractal) dimensions larger than four, such as Cartesian 
products between Euclidean lattices and dual Sierpin- 
ski gaskets (for the precise definitions of these structures 
and of the fractal dimension considered here, see later 
in this article). In this way we are able to show that 
for translationally invariant structures with high dimen- 
sions, Trw,s(t) displays sharp peaks, while for fractals or 
low-dimensional structures the peaks are low and broad 
so that the quantum walk is not particularly effective in 
the sense that there exists only a low probability TTw,s{t) 
for any t. Moreover, for any structure, we evidence inter- 
ference phenomena which give rise to a non-monotonic 
time dependence of the Trw,sit)', such effects can be sig- 
nificant and must be properly taken into account when 
considering applications. 

Interestingly, the CTQW Hamiltonian used in the 
quantum search also describes, in solid-state physics, the 
dynamics of a tight-binding particle in the presence of 
static, substitutional impurities. In this context our re- 
sults show that in regular, highly-connected geometries 
(such as the high-dimensional tori) the probability of 
finding the moving particle at the impurity site is (quasi- 
) periodic in time and that the localization can be very 
effective. 

Our article is structured as follows. In Sec. |ll]we first 
review basic principles concerning Grover's search and we 
explain how it can be implemented by means of CTQWs. 
In Sec. lIIII we describe the structures used as substrate for 
the CTQW. Then, in Sec. |IV A| we present several analyt- 
ical results, later corroborated and deepened in Sec. |IVB| 
where our numerical results are shown. In Sec.|V]we focus 
on the success probability. Finally, Sec. |VI| contains our 
conclusions and discussion. In Appendix |A[ Appendix [B] 
and Appendix [C] we report the details of our analytical 
calculations. 



II. QUANTUM WALKS AND GROVER'S 
SEARCH 

Grover's search algorithm fS] is meant to solve the un- 
sorted search problem under the assumption that there 
exists a computational oracle working as a black-box 
function able to decide whether a candidate solution is 
the true solution. Hence, the oracle knows which is the 
target among the N entries. The task is to find a tar- 
get w using the fewest calls to the oracle. While the 
classical algorithm requires exhaustive searches implying 
0{N) queries, Grover's algorithm is able to find w us- 
ing 0{y/N) queries, giving rise to a quadratic speed up 

0110]. 

A short outline of the idea behind Grover's algorithm 
in the presence of a single marked target is as follows: 
First of all, one associates each of the N index integers 
with a unique orthonormal vector |a;) = |2), |A^) in 
an A^-dimensional Hilbert space. Then, one chooses as 



initial state 

1 ^ 

I-)- 7^ El-)' (1) 

which is delocalized over the entire set of states \x) with 
equal weights at every site x. This is the least biased 
initialization one can arrange, given the available infor- 
mation (since each of the N nodes is in principle equally 
likely to be the target index, the initial state is prepared 
as an equally weighted superposition of all N indices). 

Now, to perform the search, one needs to make the 
state \s) evolve into a state that has almost all the am- 
plitude concentrated in just the \w) component. In such 
a way a single final measurement will find the system 
in the state \w), hence revealing the identity of the tar- 
get index. More precisely, we need an evolution operator 
whose repeated application makes the amplitude of \w) 
grow with the number of iterations [6]. This task was 
originally accomplished within the standard paradigm 
for quantum computation [8^, namely using a discrete 
sequence of unitary logic gates, while in the last years 
several different implementations have been introduced 
in which the state of the quantum register evolves con- 
tinuously under the influence of a driving Hamiltonian 
[l3l [2n [22] . In particular, here we follow the approach 
developed by Childs and Goldstone [12] which relies on 
CTQW [IT]. As already mentioned, due to their ver- 
satility, CTQWs allow to model any discrete database. 
In fact, a generic discretizable database can be repre- 
sented by a graph Q = {V, E} made up of a set of nodes 
V = {1,2,..., N}, each corresponding to a different item, 
and of a set of links E joining nodes pairwise in such a 
way that the topology of the graph mirrors the arrange- 
ment of the database. The graph Q can be algebraically 
described by the adjacency matrix A whose entry Aij 
equals one if nodes i and j are connected, otherwise it is 
zero (also for the diagonal elements). From A one can di- 
rectly calculate the degree matrix Z, which is a diagonal 
matrix with elements Zij = ZiSij, where Zi = X]jeAf^»i 
is the coordination number (or degree) of the i-th node, 
that is, the number of its nearest neighbors. 

CTQWs on a graph are defined by the Laplacian ma- 
trix L = Z — A and obey the following Schrodinger equa- 
tion for the transition amplitude ak.j(t) from state \ j) to 
state \k) [TT]: 

d ^ 

-a^At) = ~iY,Hkiai,j{i)^ (2) 
1=1 

where the Hamiltonian is just given by H = L, therefore, 
the time is dimensionless and given in units of the cou- 
pling elements H^i. The formal solution for aujit) can 
be written as 

afe,,(t) = (fc|exp(-zm)|j), (3) 

and its magnitude squared provides the quantum me- 
chanical transition probability iTj^ jit) = |afcj(t)p. No- 
tice that L is symmetric and non-negative definite; its 
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ground state, corresponding to eigenvalue 0, is given by 
\s)- 

Due to the unitary time-evolution generated by H, 
CTQWs are symmetric under time-inversion, which pre- 
cludes TTkj(t) from attaining equipartition at long times. 
This is different from the behavior of classical continuous- 
time random walks. Moreover, CTQWs keep memory of 
the initial conditions, as exemplified by the occurrence of 
(quasi-) revivals [S^ 1^ . 

Now, the "oracle Hamiltonian" is given by 

n^ = -\w){w\, (4) 

whose ground state, with energy —1, is just the marked 
item \w); all other states have energy zero [21]. Then, 
the Hamiltonian H governing the time evolution of the 
quantum walk is 

H = 7L + H^, (5) 

where 7 is a tunable parameter with units of inverse time, 
hence dimensionless here. 

The success probability t^Z si^) defined as (see also 

m) 

7rlJt)^\{w\eM^^m\-^)\', (6) 

i.e., as the transtion probability to be at time t at the 
target w when starting in the delocalized state | s) . 

Here we study its dependence on time t and on the pa- 
rameter 7 and we especially aim to evidence the existence 
of an optimization parameter 7max, possibly depending 
on time, which maximizes 7r2,j(t). Notice that, having 
fixed 7, the time t can be measured in terms of the num- 
ber of queries of the discrete Grover oracle [IB] . 

III. FRACTAL AND HIERARCHICAL 
STRUCTURES 

Before proceeding, it is worth introducing the geo- 
metrical structures on which we are focusing, first the 
dual Sierpinski gasket (DSG), the T- fractal (TF), and 
the Cayley tree (CT); examples of these structures are 
shown in Fig. [T^, Fig. [T]3, and Fig. [T]:, respectively. No- 
tice that these structures differ significantly from those 
analyzed previously: While hypercubic lattices with pe- 
riodic boundary conditions (toroids) are translationally 
invariant, the aforementioned structures are not. 

The DSG, TF and CT can be built iteratively; at the 
g-th iteration we have the fractal of generation g (see e.g. 
[ini IMl ISO])- The DSG and the TF are examples of ex- 
actly decimable fractals for which the fractal dimension, 
df, and the spectral dimension, d, are exactly known. 
While the fractal dimension relates the number of nodes 
inside a sphere to the radius of the sphere [31] , the spec- 
tral dimension is obtained from the scaling of the eigen- 
modes of a given structure (phonons for lattices, fractons 
for fractal substrates) [32j , most simply seen in the prob- 
ability of a random walker to be (still or again) at the 
original site [55| . 



Here, we have namely df = In 3/ In 2 « 1.585 and 
J= 21n3/ln5 « 1.365 for the DSG, and df = In 3/ In 2 
and d = 2 In 3/ In 6 « 1.226 for the TF. We recall that 
for the usual, translationally invariant lattices the spec- 
tral and fractal dimensions coincide with the Euclidean 
dimension d, namely d = df = d = \ for the chain, 
d = df = d = 2 for the square lattice and so on. On the 
other hand, on fractal structures d often replaces d when 
dealing with dynamical and thermodynamical properties 
[?7] . Also, for fractals d < df and df is smaller than 
the Euclidean dimension of the embedding space, that 
is, 2 for DSG and TF. The CT is no fractal in the clas- 
sical sense, since its growth with increasing generation is 
exponential. However, the CT is built in a hierarchical 
manner, analogous to the TF. We also notice that both, 
the TF and the CT, are trees (and are hence devoid of 
loops) and exhibit a large number 0{N/2) of end nodes. 

As we will show in the following, CTQW on the DSG, 
the TF and the CT can display a low probability 7r2,^s(i) 
for any t and 7 when compared with the case of the 
translationally invariant lattices, especially when the di- 
mension of the lattice becomes larger than 4. Therefore 
one can ask whether the probabilities ttJ, ^ (i) can be im- 
proved by adopting (fractal or hierarchical) substrates 
which display a large spectral dimension and a large co- 
ordination number. 

For instance, we can build up such structures by com- 
bining the DSG and Euclidean lattices by means of 
Cartesian products. In general, the Cartesian product 
of two graphs Gi = {Vi,Ei} and G2 = {^2,-E'2} is a 
graph G = Gi >^ G2 with the vertex set Vi XV2, and such 
that two nodes {xiX2,yiy2) are adjacent if {xi,yi) G Ei 
and X2 = y2, or xi = yi and {x2,y2) S ^'2- It has been 
shown |35j that the spectral dimension d on the product 
graph G is then the sum of the corresponding dimensions 
of the two initial graphs di and d2- 

We combine in this way the DSG with the chain Li, 
with the square lattice L2 and with the cubic lattice L3, 
to obtain more complex structures displaying spectral di- 
mensions approximately equal to 2.365, 3.365 and 4.365, 
respectively. 

Finally, it should be underlined that, dealing with such 
structures, the location of the target, i.e. the node lu, 
also (quantitatively) matters. In the numerical analy- 
sis of Sec. |IVB| we place the target on a peripheral site, 
which means, without loss of generality, the apex for the 
DSG, the leftmost site for the TF and an outmost site for 
the CT (see Figjl]). We expect that a peripheral position 
for the target site does not correspond to an optimal sit- 
uation and, since a priori the target position is unknown, 
this choice prevents an overestimation of the probability 
of finding the target. 
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a) 




FIG. 1: Examples of fractal structures considered here; in 
general the star indicates the position of the trap. Panel a: 
Dual Sierpinski gasket of generation 4 and volume A'^ = 3^ = 
81; due to symmetry, the three corners are equivalent. Panel 
b: T-fractal of generation 3 and volume A'^ = 3^ + 1 = 28; 
due to the symmetry all outmost sites (at distance 2®~^ from 
the central node) are equivalent. Panel c: Cay ley tree of 
generation 5 and volume N = Sx2^ — 2 = 94; all of the 
N/2 + 1 outmost sites are equivalent. 



IV. OVERLAPS AND TRANSITION 

Let us consider the Hamiltonian of Eq. ([s]) and denote 
the corresponding set of eigenstates and eigenvalues by 
{IV'fc)} and {Ek}, respectively. Now, for large 7 the con- 
tribution of Hu, to H is negligible and the ground state 
Itpo) is close to \s). On the other hand, as 7 — > the 
ground state is close to \w) since the weight of L is small 
and, from perturbation theory, we expect that |s) is close 
to lipi), i-e. to the first excited state of H [T^ . 

As pointed out by Childs and Goldstone [T^, there 



exists an intermediate range of 7 where, for complete 
graphs and hypercubic lattices with dimension d > 4, \s) 
switches from the ground state IV'o) to the first excited 
state IV'i); in the very same region of 7 the target state 
\w) switches from a state with large overlap with to 
the ground state IV'o)- Therefore, by varying 7, we can 
find a particular value for 7 for which the Hamiltonian 
H can evolve the state \s) into a state close to the target 
state \w). The existence and the narrowness of such a 
range for 7 are crucial for a large success probability. 

A. Analytical results 

As anticipated, the CTQW under study can yield ef- 
fective results if the Hamiltonian H is able to rotate the 
state \s) to a state with a large overlap with \w). For 
this to occur a first condition which needs to be fulfilled 
is that there exists an intermediate range of 7 over which 
the state |s) has a substantial overlap with both {ipo) and 
In this subsection the occurrence of such a condi- 
tion is investigated analytically for an arbitrary struc- 
ture, following arguments similar to those exploited in 
[T^ for d-dimensional hypercubic lattices; the details of 
the calculations are given in Appendixes A and B while 
here we just report the main results. 

First of all, we define \(f>k) and £{k), the fc-th eigenstate 
and eigenvalue of the Laplacian L, respectively. In the 
basis of the eigenstates the target site |w) can be 
written as 

N-l 

\w) = ^ ak\(pk), (7) 

fc=0 

where ak = {w\4>k)- 

As derived in Appendix |Aj once S,j = 
J2k^Q hfeP/[^(^)]"' is set, the overlap of \s) with 
the ground state jV'o) or with the first non-degenerate 
excited state jV'i) is (depending on 7) limited from below 
by the same bound. We find namely 

l>l(.|^o)P>l-^^(^, (8) 
for 7 > <^i, and also 

l>l(^l^i)P>l-^^(^f^, (9) 

for 7 < ^1. 

Now, we define the (size dependent) critical value 7 as 
that value of 7 for which 

\{s\^o)\' = \{s\^i)\' (10) 

is satisfied, i.e., the 7- value for which the projection of 
state s onto the ground and the first excited state has the 
same magnitude. This is particularly interesting in the 
limit TV — 00, when a level crossing from state |'0o) to 
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state can occur. In fact, we find that if C2/[^(Ci — 
7)^] converges (at least in the Umit TV — >■ 00) to as 
7 approaches ^1 (from different sides), then KsjV'o)! f^nd 
|(s|?/'i)| both approach 1 (see Eqs. ([s]) and namely 
a transition occurs at 7 = 1^1. Notice, however, that the 
condition for this to occur is in general non-trivial as 
^1 and ^2 both depend on N. In Appendix |X] we find 
a sufficient condition in the Laplacian spectrum and, in 
particular, in Appendix[B]we prove that such a condition 
holds for the DSG for which £{k) is exactly known [26l 
[251 ; in this case we find that 

6 < CN°'+^/^, (11) 

and 

6 < CAf"+^/^ (12) 

where a is a parameter depending on the given net- 
work. In general, — 1 < a < 0; for hypercubic lattices 
a = —1, regardless of their dimension, while for the DSG 
we can numerically estimate a as being a ~ —0.9 (see 
Appendix [A|) . Therefore, we expect (for DSG) 7 to be, 
approximately: 

7 « CiV"+2/'^~. (13) 

This result is consistent up to logarithmic corrections 
with the critical points found in [12 for the linear chain 
and the square lattice, namely 7 ~ iV and 7 ~ logiV, 
respectively. 

Finally, we point out that the critical parameter 7 pro- 
vides interesting information in the context of quantum 
adiabatic computation [13l |34j : 7 represents a threshold 
below which we can expect \w) to have a large overlap 
with the ground state. 

B. Numerical results 




FIG. 2: (Color online) Overlaps for a DSG of generation g = 3 
(up) and g — 6 (bottom) as a function of (the dimensionless) 
7, see text for details. 



We now consider the three examples of low- 
dimensional inhomogeneous structures described previ- 
ously, for which the overlaps of the initial state \s) and 
of the target state \w) with and jV'o) SLie shown in 
Figs. [2]|4] for different generations g, as a function of the 
parameter 7. These plots evidence that there exists an 
intermediate range of 7 where the state \s) changes from 
having a large overlap with the first excited state to hav- 
ing a large overlap with the ground state. In the same 
region of 7 the overlap is significant for struc- 

tures of small size (top panels in Figs. [2]|4]), while it is 
still very small when the size is large (bottom panels in 
Figs. 2]|4|. This is vastly distinct from the situation found 
for hypercubic lattices IT5t where close to 7 the overlap 
Kwli/ji)!^ is significant. 

By comparing the plots obtained for the DSG (Fig. [2]), 
the TF (Fig. [3]), and the CT (Fig. we notice that 7 
depends sensitively on the underlying topology. In fact, 
going from DSG to CT and then to TF we notice an 



amplification of the transition region, which is for TF 
most spread out, requiring relatively large values of 7 in 
order for |s) to have a large overlap with the ground state. 
Such effects can be ascribed to the absence of loops and 
to the existence for TF of a large number of peripheral 
sites scattered throughout the whole structure which give 
rise to localization effects [H]. 



We now calculate 7 according to Eq. ( 10 1 and for sev- 



eral values of g; numerical data and relative best fits are 
depicted in Fig. [5] For the DSG, numerical points are 
best fitted by the function 7 « 30 559 ^0.55^ ygj.y 
good agreement with the analytical findings. In fact, ac- 
cording to our analytical investigations, 7 sc ales with the 
size of the database like A^^+^Z^^ (see Eq. (13l), where, 
for the DSG, a + 2/d ^ a + log 5/ log 3 « 0.57 (a is 
taken to equal —0.9). Let us now consider the case of 
the TF: From numerical data the best fit turns out to 




FIG. 3: (Color online) Overlaps for a T-fractal of generation 
g = 3 (up) and g = 6 (bottom) as a function of 7. The 
symbols are as in Fig. [2] 



FIG. 4: (Color online) Overlaps for a CT of generation g = 3 
(up) and g = 10 (bottom) as a function of 7. The symbols 
are as in Fig. [2] 



be 7 w 3° ™f ArO-70_ interestingly, this result is still 
in very good agreement with the analytical approxima- 
tion of Eq. (13). In fact, for the TF the exponent gets 
a + 2/d = a + log6/log3 ~ 0.73 (where, again, a is 
taken to equal —0.9, consistently with the estimates given 
in Appendix [a]). Such a consistency might suggest that 
Eq. ([T3| is valid not only for the DSG but for any (exactly 
decimable) fractal with d < 2. As for the CT, not being 
a fractal, Eq. ( 13 ) does not hold. Indeed, we find that 



the value of 7 corresponding to Eq. ( 10 ) grows linearly 



with the generation of the fractal, namely logarithmically 
with N. In Fig. [5] numerical data are fitted by the curve 
J — g — 1 (notice the semilogarithmic plot) . 

Apart from this, the plots shown in Figs. 2]|4 look 
rather similar. In particular, for networks of large sizes 
|(w|'0o)P decays with 7 more rapidly than |(s|'(/'i)p. 
Hence, the range of 7 over which the transition occurs is 



wide, analogously to what happens on low-dimensional 
hypercubic lattices (see [H]). As shown in the next sec- 
tion, this has important effects on the behaviour of the 
success probability and suggests that a sharp transition 
is associated with an effective search algorithm. 

In order to sketch the role of the position of a target 
we show in the inset of Fig. [4] for the CT of (7 = 10 
the case of a target placed on a nearest neighbor of the 
central node. We see that the transition region is shifted 
towards lower values of 7. This means that a more central 
placement of the target is improving the probability for 
the CTQW to reach the target. Analogous results were 
also found for DSG and for TF. Thus, limiting our focus 
to peripheral nodes will prevent us from overestimating 
the success probabilities. 

We now consider fractal structures exhibiting large 
spectral dimension; in particular, we focus on fractals 
obtained from Cartesian products, such as DSG x Li, 
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FIG. 5: (Color online) Estimate of 7 for DSG (•), TF (*) 
and CT (o) in a semi- logarithmic scale. The continuous lines 
represent the best fits. 



DSG X L2 and DSG x L3, as introduced in the previ- 
ous section. Again, we place the target on a "peripheral 
site" , i.e. on a minimally connected site; this displays the 
coordination number 5, 7 and 9, for DSG x Li, DSG x 
L2 and DSG x L3, respectively. We calculate for these 
structures the overlaps of the initial state \s) and of the 
target state |w) with ji/'o) and IV'i); results for DSG x 
L2 (J = 3.365) and for DSG x L3 {d ^ 4.365) are shown 
in Fig. 6. 

As stressed at the beginning of this section, these plots 
provide some information about the sharpness of the 
transition from state |s) to the state \ipo) and from state 
\w) to the state IV'o)- However, around 7 also a signifi- 
cantly large overlap |(?i'|V'i)P is required. 

Here, the transitions are still rather smooth although, 
by increasing d, the region of 7 values, over which the 
curves representing [(u'li/'o)!^, KsIi/iq)!^ and |(s|V'i)P in- 
tersect, is shrinking. On the other hand, the overlaps 
between \w) and the first excited state \tpi) are negligible 
for all values of 7. 



V. SUCCESS PROBABILITY 

We now turn to the success probability 7r2,_s(i), Eq. (j6]), 
and we investigate numerically its dependence on t and 
on 7. Because of its dependence on time, tt^, g{t) carries 
more information than the previously discussed overlaps. 

We first analyze the case of complete graphs and of d- 
dimensional hypercubic lattices (for which the time de- 
pendences of 7r2, 5 (t) have already been determined for 
special choices of 7 in Ref. [12]) before turning to the 
DSG, the TF and the CT. 

We start our analysis from the complete graph Kj^ for 
which, as shown in [T^], at 7 = 1/N the ground state 
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FIG. 6: (Color online) Upper panel: Overlaps for a DSG 
x L2 structure given by the Cartesian product of a DSG of 
generation g — 4 and a square lattice of size L = 8; Bottom 
panel: Overlaps for a DSG x L3 obtained from a DSG of 
generation g — 4 and a cubic lattice of size L = 4. The lines 
are as in Fig. [2] 



changes sharply from \s) to \w). This transition takes 
place a,tt — tt^/N /2. Due to the special topology of if at, 
we are able to calculate 7r2, g{t) exactly, obtaining 
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4{N - 1) sin^ {t{y/4'y + (N-f ~ 1)^/2) 
4 + 7(7V- 1/7)2 

(14) 

see Appendix [C| for details. 

In Fig. [7] we show '!Tw,sit,j) for the complete graph 
with N = 124 and N = 3125. We evaluated the figures 
both numerically, by first diagonalizing H in Eq. ^ and 
projecting on the states |w) and \s) and also by mak- 
ing use of Eq. ( 14 ) . The results are numerically indis- 
tinguishable. From Fig. [7] we see that, around the val- 
ues 7 = 8 x 10-^ and 7 = 3.2 X 10"^ for TV = 124 



200 400 600 800 

t 



50 100 150 200 

t 



FIG. 7; (Color online) Contour plot of the success probability 
TTw,s{t,j) as a function of (the dimensionless) time t and of 
7 for the complete graph of size N = 124 (top) and A'^ = 
3125 (bottom). One can notice that for 7A'^ = 1, namely 



7 = 8.1 • 10 
a period r = 



(top) and 7 = 3.2 ■ 
rv/i^ « 35 and r ? 



10 (bottom), TrZ,s{'t) has 
i 176, respectively. 



and N = 3125 respectively, TrZ.sii) reaches values very 
close to 1. In fact, analyzing Eq. (14) (see Appendix [C|) , 
one finds that tt^ g (t) attains its maximal value of 1 for 
7iV = 1 and for t = Try/N{k + 1/2), where k e Z. There- 
fore, 7max = l/N. Furthermore, due to the fact that the 



period between maxima is T = n\/ N for 7max, it follows 



that the CTQW takes 0{vN) queries to find the target, 
in agreement with previous results [121 [21] . On the other 
hand, the exact dependence on t and on 7 also allows to 
highlight the oscillating behaviour of ttJ ^ (i). This means 
that, although we properly select 7max ~ '^/N, the result 
for the walk depends sensitively also on t. In particular, 
for t = kn-'/N, k G Z, the success probability is minimal 
and equals l/N (which also corresponds to the absolute 



FIG. 8: (Color online) 5-dimensional torus with linear size 
L = 5; Top: Overlaps (symbols are as in Fig. [2|; Bottom: 
Contour plot of the success probability n^^sit, 7) as a function 
of t and of 7; the dashed white line represents 7 « 0.12. 



minimum) . 

For the hypercubic lattices the overlaps of the states \s) 
and with and |-0i) have already been analyzed 
in [12]; here we display analogous plots (but for larger 
sizes) in order to compare them with the corresponding 
success probability tt^, g{t). In Fig.lslwe consider the case 
of a 5-dimensional torus (i.e. a nve-dimensional cubic 
lattice with periodic boundary conditions) of linear size 
L = N^/'^ = 5: The transition at 7 « 0.12 is very clear 
(see the top panel) and the success probability is sharply 
peaked just at 7max ~ 7 (see the bottom panel). 

Away of the critical point 7 the success probability 
quickly decays as a function of 7: it is just in the region 
of largest overlap between the initial state |s) and the 
target state \w) (namely around 7) that one expects an 
optimal success probability. 

Again, we notice that 7r2,s(<) oscillates in time; for 
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FIG. 9: (Color online) 2- dimensional square torus of linear 
size L = 56. Top: Overlaps (symbols are as in Fig.[2]|; Bot- 
tom: Success probability as a function of time and 7; the 
dashed white line represents 7 « 0.67. 



FIG. 10: (Color online) Contour plot of 7rJ,^(t) for the DSG 
of generation 3 = 3 (top) and 3 = 6 (bottom); the "critical" 
values 7 are represented by the dashed line: 7(3) ~ 1.20 and 
7(6) ^ 7.38 (see Fig.|2|. 



L — b and 7 — 7max the success probability ranges from 
about to about 0.8. Moreover, for a given time t, 7r2,_s(t) 
decays very fast as I7 — 7max| increases. For instance, 
TfO ■''7— (70)/7rT;»-(70) « 0.05. As a result, the computa- 
tional procedure for this structure can be very efficient, 
provided that the parameter 7 can be sensitively con- 
trolled. 

For hypercubic lattices of dimension d = 4, d = 3 and 
d — 2 (only the latter case is depicted in Fig. |9]) the 
peaks get more and more broadened and are of smaller 
magnitude. Notice that the low peaks obtained are in 
agreement with the analytical results found in |12j which 
predict for the 2-dimensional torus a vanishing success 
probability for N ^ 00. 

We now turn to the analysis of the success probability 
for the DSG, the TF and the CT, represented in Fig- 



ures 10 11 and 12, respectively; the upper panels display 



the situation for small, the lower panels for larger net- 
works. First of all, we notice that, as previously found 
for regular lattices [H], also for DGS, TF and CT, nl_^{t) 
exhibits peaks which are lower and lower as the size N 
is enlarged. Moreover, for small sizes (upper panels) 
7max ~ 7, namely the maxima for the success probabil- 
ity occur for values of 7 which are approximately equal 
to 7. On the other hand, for large sizes (lower pan- 
els), in the temporal range considered here, 7 provides 
an upper bound for 7max- However, the most striking 
feature which emerges from the comparison between the 
contour plot of the success probability for translation- 
ally invariant structures (see Figs, [s] a nd [9| ) and for frac- 
tal/hierarchical structures (see Figs. |10|-12) is that for 
the latter n^^sit) is much more rough and broadened. 

Now, it is worth comparing the 5-dimensional torus of 
Fig. |8) the square torus of Fig. |9] and the CT of Fig. 12 
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FIG. 11: (Color online) Contour plot of ttZA^) ^he TF 
of generation g = 3 (top) and g ~ 6 (bottom); notice that 
7(3) « 2.59 and 7(6) « 21.31 (see Fig.|3|. 




20 40 60 

t 




FIG. 12: (Color online) Contour plot of irZ^sit) for the CT 
of generation <? = 3 (top) and g — 10 (bottom); notice that 
7(3) « 2.06 and 7(10) ^ 8.87 (see Fig.|4|. 



since all three are of comparable size N. From the com- 
putational point of view the 5-dimensional torus corre- 
sponds to the best situation: nj, ^ {t) is sharp and reaches 
its maximum value around 0.8 after approximately 100 
unit steps; the CT displays a success probability around 
0.12 after 150 time steps; the 2-dimensional torus cor- 
responds to an ineffective candidate situation: after 300 
time steps the peak is still less than 0.1. 

In order to compare also to structures with spectral di- 
mensions larger than four, we show in Fig. 13 the success 
probabilities tt^, ^(i) for the cartesian products of a DSG 
of generation g = 3 and a square lattice of size L = 8 as 
well as a cubic lattice of size L ~ 8. Although the max- 
ima of the success probabilities are in both cases larger 
than the ones for the structures with low dimensions, 
they show still a fairly unregular pattern. This is in con- 
trast to the highly regular structure of the 5-dimensional 



torus, see Fig. |8] 



Finally, we stress that, as a result of interference phe- 
nomena, TT^, g{t) oscillates with time. This has some im- 
portant consequences: Although we can determine and 
set the optimal 7max, the probability of the CTQW reach- 
ing the target depends on the instant of time at which 
it is calculated. In particular, oscillations are "faster" 
for systems of smaller size; for example for the complete 
graph we find a period r(7) = 27r/-\/ {jN — 1)^ -I- 47 (see 
Eq. (14|), namely T(7max) = tt\/N, while for the 5- 



dimensional torus the numerical analysis makes it pos- 
sible to estimate a period r which grows exponentially 
with the lattice size L; for L = 5 we get r w 200 (see 
Fig.|8|. 
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FIG. 13: (Color online) Contour plot of TrZ^sit) for the carte- 
sian product of a DSC of generation g = 3 and a square lattice 
of size L = 8 (top) and a cubic lattice of size L — 8 (bottom). 



VI. CONCLUSIONS AND PERSPECTIVES 

In this work we considered CTQWs mimicking 
Grover's quantum search problem; we especially focused 
on how the topology of the space over which the walk 
takes place affects the position and sharpness of the tran- 
sition of the ground state. Previous studies [T^ high- 
lighted that for translationally invariant structures, such 
as the hypercubic lattices, the quantum walk can be 
highly efficient for sufficiently high dimensions, i.e. d > A. 
However, here we evidence that on generic graphs the di- 
mension does not represent the key geometric parameter; 
indeed, both the (average) coordination number and the 
fact that the structure is translationally invariant or not 
determine the sharpness of the transition. In fact, on 
the one hand, a high coordination reduces the distance 
among nodes and increases the possibility of interference 



effects, on the other hand, (in the absence of a target) 
translational invariance prevents the emergence of local- 
ization effects [29] . 

In particular, we considered the success probability 
7r2,^(t) (here |s) and \w) are the initial and the tar- 
get state, respectively) as a function of the computa- 
tion time t and of a properly tunable parameter 7 en- 
tering the Hamiltonian. We showed that for highly di- 
mensional (d > 4) translationally invariant structures 
(5-dimensional torus) there exists a narrow range of 7 
around a specific value 7 where the ground state |'0o) un- 
dergoes a transition from having a large overlap with \s) 
to a state with a large overlap with \w). This corresponds 
to a sharply peaked success probabihty: n^j^sit) displays 
a set of maxima just at 7; the first one is reached after 
0{-\/N) queries. Conversely, for structures with low coor- 
dination number and/or fractal or hierarchical topology - 
such as the cubic (L3) and square (L2) tori, the DSG, the 
TF, the CT and the Cartesian products DSG x L2 and 
DSG X £3 - the transition from the initial state to the 
target state takes place over a wider region of 7 around 
the value 7. As a consequence of such a spread-out tran- 
sition, the success probability displays broadened peaks 
whose locations depend on time. Therefore, for the non- 
translationally invariant structures considered here, even 
with large fractal dimension [d > 4), the large success 
probabilities found for high-dimensional periodic lattices 
are not recovered. These results, in agreement with pre- 
vious findings pR|, highlight a possible connection be- 
tween the sharpness of the transition occurring at 7 and 
the efficiency of the search algorithm. A mathematical, 
rigorous proof stating whether a sharp transition is a nec- 
essary condition for a good algorithm, which is beyond 
the aim of this article, could provide a very useful tool 
for further investigations on quantum search algorithms. 

For the DSG we also proved that 7 scales like iV^/''+" 
(— 1 < a < 0); interestingly our results suggest that such 
a scaling might be generalized to all (exactly decimable) 
fractals with spectral dimension d < 2. 

Apart from the deterministic fractals considered here, 
it will be extremely interesting to also consider disor- 
dered structures such as percolation clusters and ran- 
dom graphs characterized by a degree distribution P{z). 
These networks display a tunable average degree (z), 
which, in principle, can take values ranging from for 
totally disconnected networks up to N for completely 
connected networks. According to the results discussed 
here, we expect that for a sufficiently large {z) and for a 
sufficiently peaked P{z) the transition from the ground 
state \s) to a state with significant overlap with \w) oc- 
curs sharply and that, consequently, the CTQW speeds 
up. 



Acknowledgements 

Support from the Deutsche Forschungsgemeinschaft 
(DFG) and the Fonds der Chemischen Industrie is grate- 



12 



fully acknowledged. EA thanks the Italian Foundation 
"Angelo della Riccia" for financial support. 



Appendix A: Generic structures 

Let us denote by and by {£{k)} the sets of eigen- 

states and of eigenvalues of the Laplacian L, respectively. 
On the basis of the eigenstates the state \w) localized at 
the target site can be written 



N 

E 

fc=0 



ak\4>k)- 



(Al) 



Now, let us consider the complete Hamiltonian H and the 
corresponding eigenvalue equation for the state labeled a: 



(A2) 



As shown in 12 , when one sets Ra = |(w|?/'a)P, it is 
possible to write 



10" 



X 



10" 







\ jV-0.9 


^DSG 












10^ 


10^ 


10^ 



N 



FIG. 14: Numerical estimate of a for the DSG and the TF. 



whose derivative is 



and 



(A3) 



N 



F'{E) ^ 



\ak\ 



k=l 



{^£{k)-EY 



(A9) 



F{Ea) 



1 



1^-Ea 



1. 



(A4) 



Here we notice that Eq. |A3[ and therefore also Eq. |A4[ 
hold when the eigenvalue Ea is non-degenerate (for ex- 
ample, this condition is not fulfilled if we place the trap 



on the central node of the TF and CT) . Equations ( A3 ) 
and (A4) allow then to express the overlap of a given 



state |i/'a) (with non-degenerate Ea) with \s) as 

1 Ra 



\{sHa)Y 



ElF'{Ea)N NEr 



where 



and, in particular 
F'iEa) - 



(7L - Ey 



(7L-K)^ 



w 



1 



(A5) 



(A6) 



(A7) 



(see |12| for more details). It is worth underlining that 
Eq. (A5) found in [T^] for hypercubic lattices actually 



holds for all structures. The topological details are con- 
tained in Ea and F'{Ea)- 

We now proceed with the calculations without making 
any assumptions on the topology of the database. First, 
by using Eq. (m) and OfcO^ — la^p, we rewrite F{E) as 



HE) = E 

k 



akttk 



-i£{k)-E 



E 

k 



-f£{k)-E' 



(A8) 



For a generic index /c, < la^p < 1; while the lower 
bound is clear, the upper bound derives from the fact 
that the eigenstate can not correspond to any Lapla- 
cian eigenstate. For Euclidean lattices Bloch's theorem 



makes it possible to write {x\(f>k) = e^^^ js/N so that 
IflfcP = for any fc- On the other hand, for a 

generic connected structure, a priori, one can only write 
koP = as a consequence of the fact that the Lapla- 
cian eigenstate corresponding to the smallest eigenvalue 
f (0) =0 is just |0o) = |s)- This suggests a proper restric- 
tion of the previous upper bound: max^^o |a/cP ^ 
with C > and — 1 < a < 0, both depending on the 
particular topology chosen. In particular, for Euclidean 
structures, C = 1 and a = —\ and one expects that 
the more inhomogeneous the topology, the larger a. By 
means of numerical calculations we can estimate a: for 
not too small DSG and TF we find that a ~ —0.9 (see 
Fig. 14). 

Now, before going on, we define the quantity ^j, which 
will be useful in the following: 



1 



(AlO) 



where we set X^fc^^^o ^ /['^(^)P = ^^"^ used the upper 
bound Iflfcp < CA^" < 1; for hypercubic lattices = 
C^jjN which can be approximated by an integral [12j . 



Using Eqs. (A5| and (A9) the overlap of |s) with the 
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ground state turns out to be 



{j£ik) + \Eo\r 

2 



fe^O 



{i£ik) + \Eo\y 



> 1 



(All) 



where in the first inequafity we used that the sum is pos- 
itive, while in second inequafity we used that both £{k) 
and \Eq\ are positive. From Eq. (All) and Eq. (AlOl, 
we have 



IWo>|^>l-^ 



6 



(A12) 



In general, as 7 is varied, Eq is bounded as < |£'o| < 1 
and the bounds can be improved by exploiting the fol- 
lowing 



1 = FiEo) 



\ao\ 



E, 



01 



E 



\ak\ 



l£{k) + l^ol 



\ak\ 



1 



N\E^\ ^^„ 7f(fc) N\Eo\ ' 7 



In fact, we get 



(A13) 



(A14) 



where for the lower bound we used that the first sum 



appearing in Eq. (A13| is positive and that oq = Oq = 



Now, from Eq. ( A12) and Eq. ( A14| it follows straight- 
forwardly that 



l-|(«IV'o)P< 



6 



< 



(A15) 



Following analogous arguments we find for the first non- 
degenerate excited state labelled as 1, being Ei non de- 
generate. 



\{s\i,,)\'>l-NElY, 



\ak\ 



k^ 



{^£{k)-E,f 



and 



El < 



A^(a-7)' 



from which we get 



l-\{s\^Pi)\'< 



6 



CN°'-\2 



< 



(A16) 



(A17) 



(A18) 



Notice that in this case 7 < ^1 due to the fact that Ei > 
0. 



Now, by comparing Eq. (A15) and Eq. (AI8I we can 



evidence the existence of a critical value 7 such that when 
7 approaches 7 the ground state switches from \ipQ) to 
IV'i), at least in the limit — > 00. In fact, if we take 



7 = CAT^Ci, then in Eq. ( |A15[ ) and Eq. ( |A18[ ) we can 
set 'f — {C + £)N°'(i, with e > and e < respectively, 
obtaining 



i-Ksiv-o,!)!' < 



C C2 



(A19) 



Therefore, recalling that a— 1, the condition (2 — C'(Ci ), 
as A^ — > 00 is sufficient for Ksl-^o,!)!^ 1- Otherwise 
stated, as 7 approaches 7 from above (e > 0) or from 
below (e < 0), the overlap of the ground state with 
and with respectively, is close to 1. 



Appendix B: The Dual Sierpinski gasket 

The Laplacian spectrum £{k) for the DSG is exactly 
known [5^1 [251 US] and we can therefore calculate exactly 
the quantities Ci and (2, obtaining estimates for the crit- 
ical parameter 7. 

At a given generation g, the spectrum includes the 
eigenvalue 3 with degeneracy m{3,g) = (3^^^ +3)/2, the 
eigenvalue 5 with degeneracy 171(5, g) = (3^^^ — l)/2 and 
the eigenvalues A"^ and A~ stemming from the eigenvalue 
A, belonging to previous generation and both carrying 
degeneracy m{X,g — 1). For each eigenvalue A the next- 
generation eigenvalues A^ are defined according to 



, _ 5 ± V25 - 4A 

^ — 



(Bl) 



Now, it follows directly from Eq. (Bl) that any couple 
1/A+, 1/A^ sum up as 



1 1 _ 5 
A+ ^ A^ ^ A' 



(B2) 



and applying this result iteratively to all the couples mak- 
ing up the spectrum we get 

^1 ^ E^ 

^ 5J /39-J-i + 3\ ^ 5^ - 1 

= Ey(^^j+E^ 

j=0 ^ ^ j=0 

= l(-3-4x 3^ + 7x5^). 



(B3) 



As for (^2 we can implement a similar iterative procedure, 
by noticing that 



1 



1 



(A+)2 (A-)2 



25-2A 
A2 ■ 



(B4) 
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Therefore C2 is made up of two terms stemming from the 
first and second order contributions, which are, respec- 
tively. 



5^-1 (5^' - 1) + 3 

~3 2 2 

5^-1 (5^' - 1) - 1 



and 



9-1 

E 

3=0 

9-1 
E 

55 + 80 X 39 - 154 X 5f + 19 X 5^^ 
6600 



^ ^ 25^ 39-i-i + 3 ^ 25^ 39-i-i _ 1 



(B5) 



E 

3=0 



E 

-121 - 68 X 39 + 189 X 529 
19800 ■ 



(B6) 



Subtracting Eq. ( B6 ) and Eq. ( B3 1 we finally get 
1 



k^O 



[£(fc)] 



900 



(-13-14x3» + 21x59 + 6x25»). 

(B7) 

Now, recalling that N — 3^ and that the spectral di- 
mension of the DSG is d = 2 In 3/ In 5, we can write 
59 — obtain expressions for Ci and C2 as a func- 

tion of just the volume N and the spectral dimension d 
of the substrate: 



^^"3^( 



7 



, -3 - 4iV + 7iV2/'^ ) TV^/'i (B8) 

30 V y 30 ' ^ ^ 



C2 



1 

900 



-13 - 14iV -I- 21N^/'^ + 67V''/'* 



150 



(B9) 

where the asymptotic expressions hold for large N. 

We not ice th at a s N ^ 00, (^i and ^2 appearing in 
Eqs. (|B8| and (|B9| satisfy the condition (2 = C'(Ci) 
found in Appendix |A| More precisely, C2 ~ Ci and this is 



and 1— \ {s\ip 



sufficient for 1 - 

zero as A^^". In particular, from Eq. (|B8|) and Eqs. (|A15|) 



to conve rge to 



and (AI8), the state |s) is expected to switch from the 



ground to the first excited state at 
7«CAr2/d+a_ 



(BIO) 



As explained in Sec. |IV A[ the expressions found here 
for ^1 and (2 are consistent with ^1 = ^i;^/[£{k)y ^ 
]\[2j/d obtained in [T2] for lattices. Therefore, and as 
suggested by the numeric al d at a di scussed in Sec. |IVB[ 
it is plausible that Eqs. (B8), (B9) and, above all, the 
expression for the critical parameter 7 in Eq. ( |B10 ) can 
be extended to arbitrary structures of spectral dimension 
d<2. 



Appendix C: The complete graph 

Let us start from the definition of success probability 
given in Eq. (|6|. Now, the propagator U = exp(— iHt), 
can be calculated as 



1=1 



and the success probability can be rewritten as 

2 



W.S 



1 

N 



N 

E 

k=l 



Uk,v 



(CI) 



(C2) 



Hence, in order to calculate the success probability 
TrZ^sit) need to find the elements corresponding to the 
w-th column of the l-th power of the Hamiltonian. For 
complete graphs, H' displays a high degree of symmetry 
which allows its exact calculation (see for example [21] for 
a similar calculation where the Hamiltonian is provided 
by the adjacency matrix). Without loss of generality we 
can fix w = A^ so that the Hamiltonian H is 



/N-l -1 
-1 A^-l 



V 



-1 



-1 



-1 



-1 



-1 
-1 

N ~l 



-1 



A^ - 1 - 1/7 / 



(C3) 

and it is easy to see that, regardless of I . (H '),;jv = 
(H')j,Ar, with i,j < N. Therefore, Eq. (C2]) can be 
rewritten as 

<Jt) = ]^\{N~l)U,,N + UN,N\' ■ (C4) 

Now, our task is to calculate Ui,n and Un.n for which we 
need the entries (1,A^) and {N,N) of the l-th power of 
the Hamiltonian; for better readability we set uinJI) = 
(H')i,7v and aNN{l) = (H^Ar^Ar. Thus, from Eq. ([CS]) we 
derive the following recursive relations: 

aiAr(/ -f 1) = aiAr(0 - ajVAf(Oi (C5) 

and 

aNN{l + 1) = -(A^ - l)aiw(0 + {N-1 
From their combination we get 

aiNil + 2)^iN~ 7"')aiiv(/ + 1) + T^'amCO- (C7) 
whose solution is 

aiN{l)^^^[iA + By-{A-By], (C8) 



7 ^)aNN{l)- 
(C6) 



with A = A^7 - 1 and B = ^/A^~^\^. From Eq. ( |C5[ ) 
and Eq. (C8), aNN{l) is also explicitly defined: 



dNNil) 



2l 

2' B 



[{A + bY - {A-B)'] 



(C9) 
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Now, according to Eq. (C2), we can calculate Ui,n as 



1=0 

7 



= -J- {cxp [-it{A + B)/2]- exp [~it{A - B)/2]} 
B 



T tB 
= — exp-ityl/2(-2isin — ] 
B 2 

Analogous calculations lead to 
exp [-it/2A) 



(CIO) 



N 



B 

tB , tB tB\^ , 

X \ —2i^sv[Y— + lAsin— _Bcos — (pil) 



Hence, by inserting Eqs. (ClOl and (Cll) into 



Eq. ( C4 1 , and performing the summation 

47(iV- 1) 



1 



<,sii) iV + 1 + 47 - 27V7 + iV272 



X sin^ {t{y/\ + 47 - 2iV7 + 7V2^2)/2){pi2) 



which provides the exact success probability for the 
complete graph of size as a function of 7 and of 
t. We can notice that tt^, ^ {t) is strictly larger than 
zero and that it can be equal to 1 when both t ~ 
7r(2fc + l)/v^l + 47 - 2iV7 + 72iV2, with /c e Z and 
47/(1 + 47 - 2A^7 + iV2^2) ^ I g^^g satisfied. The lat- 
ter condition holds for 7A^ ~ 1, just consistent with [12j : 
there it is found that for 7iV = 1 the walk rotates the 
state from \s) to \w) and that the gap between the cor- 
responding energies is smallest. For 77V = 1 the first 
condition reads t = y/NT^ik + 1/2). Hence, for the exact 
7 and at the right time the success probability is unitary, 
but for larger volumes the right times get sparser . We 
also notice that when N is large and Nj ~ 1, Eq. (C12) 
can be simply rewritten as 



(C13) 



7TZ^,{t)^sm' (t/VN) 



From Eq. (C12 ) we deduce that for a given time i, tt^, s{t) 
decreases asJjN — 1 \ gets larger. 
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